Ferromagnetic transition temperature enhancement in (Ga,Mn)As semiconductor by 

carbon co-doping 
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We present a theoretical study of (Ga,Mn)(As,C) diluted magnetic semiconductors with high C 
acceptor density that combines insights from phenomenological model and microscopic approaches. 
A tight-binding coherent potential approximation is used to describe the electronic structure in the 
presence of Mnc a and Ca s impurities. We find only a small effect of C on the distribution and 
coherence of electronic states close to the top of the valence band and on the coupling between 
Mn moments, even at doping levels of several per cent. These results justify applying the model 
of ferromagnetic Mn-Mn coupling mediated by itinerant holes in the valence band also to C doped 
samples. The increase of ferromagnetic transition temperature due to the presence of C acceptors is 
illustrated by calculations that use the kp Kohn-Luttinger description of the GaAs valence band and 
assume systems where Mn local moment and itinerant hole densities can be varied independently. 

PACS numbers: PACS numbers: 71.20.Nr, 75.30.Et,75.50.Pp 
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Prospects for new device functionalities in all- 
semiconductor spin-electronic structures rely on the re- 
alization of a ferromagnetic semiconductor operating at 
room temperature. An important milestone in this ma- 
terial research was the discovery five years ago 1] of fer- 
romagnetism in Mn-doped GaAs with the Curie tem- 
perature T c — 110 K [J. In the Gai-^Mn^As diluted 
magnetic semiconductor (DMS) with x m 1 — 10%, Mn 
substituting for Ga provides a local moment 5 = 5/2 
and a delocalized hole |3, |J] • Ferromagnetic coupling be- 
tween Mn moments is mediated by the itinerant holes 
via a kinetic-exchange interaction |5| . Recent progress in 
low-temperature MBE growth and post-growth anneal- 
ing techniques has led to the increase of the transition 
temperature in (Ga,Mn)As DMS's by nearly 50 K [(|. 
This success is attributed to a smaller concentration of 
carrier and moment compensating defects, especially in- 
terstitial Mn, in the optimally annealed samples Q ■ The 
latest observations are consistent with an approximately 
linear dependence of TrOTi x and hole Fermi wavevector 
predicted by theory IgL [To| . 

In this paper we address theoretically the possibil- 
ity of increasing T c in (Ga,Mn)As DMS's by a non- 
magnetic acceptor co-doping, namely by introducing sub- 
stitutional Cas impurities. Our work is partly motivated 
by a recent experimental observation of a marked en- 
hancement of the Curie temperature in Mn-implantcd 
GaAs:C samples compared to the Mn-implanted undoped 
GaAs layers After briefly discussing the differ- 

ent effects of Cas and BeQa [l3 acceptors on substitu- 
tional Mn incorporation in (Ga,Mn)As DMS's, we focus 
on two key issues related to high-density carbon dop- 
ing: (i) The effect of C impurities on the density of 



states in the semiconductor valence band is assessed us- 
ing the tight-binding/coherent-potential-approximation 
(TB/CPA) description of the disordered semiconductor, 
(ii) The strength of the hole-mediated ferromagnetic Mn- 
Mn coupling is compared for systems with and without 
C co-doping using the TB/CPA results and T c is then 
estimated from a model that combines the k • p Kohn- 
Luttinger description of the GaAs valence bands and a 
mean- field treatment of the kinetic-exchange coupling be- 
tween Mn local moments and the band- holes. 

An important issue for hole co-doping in (Ga,Mn)As 
DMS's is the change of the substitutional MnQ a forma- 
tion energy caused by the presence of additional non- 
magnetic acceptors. Although the formation energy is, 
strictly speaking, an equilibrium characteristic we as- 
sume that the basic features of its compositional de- 
pendence are reflected also in the dynamics of the non- 
equilibrium growth of the co-doped materials. A system- 
atic ab-inito study of the formation energies is beyond 
the scope of this short paper and will be discussed else- 
where |l3j . Here we present a qualitative analysis based 
on the TB/CPA calculations, described in more detail 
below. Our results suggest that Be acceptors substitut- 
ing for the same element as Mn, i.e. for Ga, lead to a 
strong enhancement of the Mnc a formation energy. At 
the same time, the formation energy of interstitial Mn 
is suppressed and, hence, an increasing fraction of Mn is 
incorporated in the form of interstitial donors or electri- 
cally neutral MnAs or Mn clusters. This scenario, which 
leads to a decrease rather than an increase of T c , has re- 
cently been established by extensive experimental studies 
of Be co-doped (Ga,Mn)As samples 12]. 

The TB/CPA calculations indicate that C acceptors 
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lead to a much weaker increase of the MriQa formation 
energy. We surmise that this property stems from the 
anomalous nature of Ca s which acts as an acceptor, yet 
has a larger Pauling's electro- negativity than As. The 
difference between the effects of BeQa and Ca s impuri- 
ties on Mnca incorporation is further enhanced at high 
doping levels where the competition of Be and Mn for 
the same lattice site starts to play a role. 

The above qualitative analysis of formation energies 
suggests that Ca s co-doping is favorable for achieving 
high Curie temperatures. Unlike typical acceptors, how- 
ever, C has a very different atomic size and, as already 
mentioned, a large electro-negativity compared to the 
atom it substitutes for. The crucial question, addressed 
in the following paragraphs, is then to what extent dop- 
ing by several per cent of Ca s changes the semiconduc- 
tor band structure and whether the model of shallow 
acceptor carrier-induced ferromagnetism still applies in 
(Ga,Mn)(As,C). 

Our calculations of the electronic structure of GaAs 
in the presence of Mnc a and Ca s impurities are done 
using the tight-binding version of the coherent potential 
approximation (see e.g. Ref. 0]). The CPA, in con- 
trast to supercell calculations, is well-suited for the mixed 
crystals with low concentrations of randomly distributed 
impurities. It provides estimates of a configurationally 
averaged density of states (DOS) and related quantities. 
The configurational averaging restores full translational 
symmetry of the lattice and makes it possible to decom- 
pose the DOS into contributions from specific points in 
the Brillouin zone. The spectral density A(k, E) then in- 
cludes a detailed information about the dispersion of the 
electronic states in the reciprocal space. In addition, the 
width of the peaks of A(k, E) defines the scattering rate 
of the band electrons in various parts of the Brillouin 
zone due to the impurities. 

The parameterization of the TB Hamiltonian provides 
a correct band gap for the pure GaAs crystal |l5j and an 
appropriate exchange splitting of the Mn d-states. Local 
changes of the crystal potential at both Mn and C impu- 
rities, represented by shifted atomic levels, are estimated 
using Ref. 0]. Long-range tails of the impurity poten- 
tials, which become less important with increasing level 
of doping, are neglected. (Note, that the Thomas-Fcrmi 
screening length is only 3-5 A for typical carrier densities 
[l7) . i.e., comparable to the lattice constant.) Also lattice 
relaxation effects are neglected within the CPA. This is 
well justified for Mnc a impurities 0,0] but becomes a 
more important issue in the case of Ca s impurities. Pre- 
vious density-functional studies found a -0.6% relative 
change of the GaAs lattice constant in a 64-atom super- 
cell with a single Ca s impurity [2(J . Although the lattice 
relaxation around Ca s may change our results quantita- 
tively, we expect that the dominant effects on the band 
structure and on the spectral broadening arise from the 
different atomic levels of C compared to As which is read- 



ily accounted for in the TB scheme. 

In Fig. H we plot the spin-polarized DOS in 
Gai_ 2; Mn a ;Asi_j,Cj, with x = y = 4% together with lo- 
cal DOS on host (As) and impurity (Cas) atoms. Despite 
the remarkable difference between the atomic levels of As 
and C, Ae p w -I eV, the local DOS on Ca s sites (thin 
full line) near the valence band edge does not differ much 
from the local DOS on the As sites (dashed line). (A still 
larger full width of the valence band than As p may partly 
explain this result.) Similarly, the total DOS (thick full 
line) is only weakly affected by the presence of C in this 
spectral region. In particular, the spin-splitting 0.26 eV 
of the valence band edge, directly related to the kinetic- 
exchange parameter J pi j, is nearly the same in systems 
with and without Ca s impurities. The shaded region in 
Fig. n shows the DOS of Mn d-states peaked near the 
energy of -4eV which is consistent with results in C free 
(Ga,Mn)As. 

It is important to point out that the substitution of C 
also has a very small effect on the line-width of A(k,E) 
for the states close to top of the valence band, i.e., elec- 
tron scattering on C impurities does not disturb sub- 
stantially the coherence of the Bloch states. This, to- 
gether with an unchanged value of J p d, implies that the 
additional disorder due to the co-doping with C should 
not have any marked effect on the carrier mediated cou- 
pling between Mn moments. We check this more explic- 
itly using the compatibility of the CPA with the Weiss 
mean-field theory. The strength of the Mn-Mn coupling 
is characterized by the energy cost of flipping one Mnc a 
moment, which can be calculated for a given chemical 
composition j^jj. In Fig. we plot this quantity as a 
function of the hole density, p, for x = y = 4% (solid 
line), and x = 4% and y = (dashed line) samples. The 
hole density is varied in the calculations indepenedently 
of x and y, i.e., we assume implicitly a compensation 
whenever p is smaller than the total density of MnQ a 
and Cas acceptors. The curves are nearly identical for 
systems with and without C, as anticipated above. This 
means that, for a given hole concentration, the exchange 
coupling is quite insensitive to the large amount of the 
additional Ca s defects. We conclude our TB/CPA con- 
siderations by stating that, as in the MnQ a acceptor case, 
the main effect of Ca s doping is a downward shift of Fermi 
energy with respect to the band edge, i.e., an increase of 
the number of holes in the valence band. In the following 
paragraphs we discuss prospects for ferromagnetic tran- 
sition temperature enhancement by adding extra holes 
into the DMS valence band. 

We use a model in which the itinerant holes are 
described by the GaAs host bands and the coupling 
to the local moments by a phenomcnological constant 
Jpd E3- This theory has been successful in de- 
scribing scmiquantitatively many non-trivial thermody- 
namic and transport properties of (Ga,Mn)As DMS's 
H ES El HI E3. In the simplest, virtual- 
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crystal mean-field version of the model, that assumes 
ferromagnetic indirect coupling between Mn ions, the en- 
ergy to flip the Mn Ga moment is proportional to the ef- 
fective field H e ff — Jpd(s). Here (s) is the mean spin- 
polarization density of the itinerant holes which increases 
with the hole density. This result is consistent with our 
TB/CPA calculations, shown in Fig. for hole densi- 
ties not too much larger than the MnQa density. Recall 
that the concentration of the substitutional Mn ions is 
given by Njtfn = 4a;/af c where ai c is the GaAs lattice 
constant (N Mn = 0.88 nm~ 3 for x = 4%, e.g.). For 
p > Nmu, RKKY oscillations of the Mn-Mn coupling |f| 
start to play a role and the increasing number of anti- 
ferromagnetically coupled Mn moments leads to a satu- 
ration or even to a suppression of T c , as seen in Fig. El 
for p > 1.3 nm~ 3 . The mean-field theory that allows 
only for collinear ferromagnetic states is therefore likely 
to break down in this high hole-density region. Note also 
that for high Nmu and low p, the direct antifcrromagnetic 
Mn-Mn interaction takes over, as suggested by negative 
Weiss exchange-field values in FigEl which sets another 
limit on the validity of the simple mean-field model. This 
low hole-density region is, however, not important for the 
high-T c (Ga,Mn)(As,C) DMS's we focus on in this paper. 

We will use now the mean-field model to estimate the 
Curie temperature as a function of the Mng a local mo- 
ment density and of the density of itinerant holes that can 
be varied independently. We emphasize that hole dop- 
ing due to C is relevant also for systems with p < Nmu 
since most of the experimental (Ga,Mn)As samples show 
some level of compensation, usually caused by intersti- 
tial Mn defects. In the calculations, the GaAs host 
band-structure is obtained from the kp, Kohn-Luttinger 
model |2^|. We neglect the ~ 10 — 20% suppression of 
T c due to spin-wave fluctuations since the effect is com- 
pensated, to a large extent, by hole-hole exchange en- 
hancement of T c [llj], also neglected in the present cal- 
culations. Details of the model are described elsewhere 
here we recall only the general T c expression 
we useMlllIS: 



knTc = 



N Mn S(S+l) J 2 pd X 
3 Omb) 



(1) 



where x is t ne band-hole magnetic susceptibility which 
is roughly proportional to p 1 / 3 . 

Fig. El shows constant-T c curves calculated for critical 
temperatures ranging from 50 K to room-temperature. 
The dashed line corresponds to p = Nmh and also indi- 
cates approximately the onset of RKKY oscillation ef- 
fects on T c . Note that the simple mean-field theory 
used to calculate the constant-T c curves is expected to 
break down above the dotted-dashed line where the large 
number of antiferromagnetically RKKY-coupled Mn mo- 
ments might lead to a decrease rather than increase of T c 
with increasing hole density. 



FigEH suggests that a substantial enhancement of the 
ferromagnetic transition temperature may be expected in 
hole co-doped samples with high Mn moment concentra- 
tion. Assuming, e.g., x — 10%, an increase of the hole 
density from p w 0.1 nm -3 to p rj 0.7 nm~ 3 leads to an 
increase of the theoretical T c from 50 K to 300 K. For 
smaller x, a larger Ap is needed to enhance T c by the 
same amount. Recent experiments have demonstrated 
that good quality (Ga,Mn)As DMS's can be grown with 
x reaching 8% |6(]. Based on the curves in Fig. [3] and our 
TB/CPA results we conclude, that (Ga,Mn)As DMS's co- 
doped with several per cent of C should not be overlooked 
among potential candidates for a room-temperature fer- 
romagnetic semiconductor. 
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FIG. 1: TB/CPA density of states in Gai-^MnxAsi-yCy 
DMS with x = y = 4%. Total DOS (thick full line), local 
DOS on host As (thin dashed line) and impurity Ca s (thin 
full line) atoms, and the DOS of Mn d-states are plotted as a 
function of energy, measured from the Fermi level. 
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FIG. 2: Energy cost of flipping one MnGa moment, obtained 
from the TB/CPA spectra, is plotted as a function of the hole 
density for x — y — 4% (solid line), and x = 4% and y = 
(dashed line) 

We acknowledge fruitful discussions with Cammy 
Abernathy, Miroslav Cukr, Vft Novak, and Dan Park. 
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FIG. 3: Mean-field, k ■ p model calculations of constant T c 
curves are plotted in the MnGa doping (x) - hole density (p) 
diagram. Dashed line corresponds to p — Nm„ = 4x/af c and 
the dotted-dashed line indicates approximately hole densities 
at which the mean-field model is expected to break down due 
to RKKY oscillations of the Mn-Mn indirect coupling. 
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